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INVESTIGATION OF ADVANCED HAZARD 
DETECTION CAPABILITIES WITH 
MULTIFREQUENCY RADAR 


Andrew L. Pazmany, Ph.D. 

ProSensing Inc., Amherst, MA 

EXECUTIVE SUMMARY 

Current state of the art Aircraft Weather Radars (AWRs) are scanning, X-band, single polarization (H-pol) and 
coherent with Doppler measurement capability. The measured radar parameters include reflectivity and Doppler 
velocity mean and standard deviation. The X-band (~9.4 GHz) operating frequency offers the optimum compromise 
between available transmit power, antenna gain with a moderate size aperture, receiver sensitivity (front-end losses 
and LNA noise figure) and acceptable attenuation in moderate liquid precipitation. The Doppler velocity standard 
deviation measurements are effective in identifying turbulence in precipitation. 

However, current X-band AWRs have the following shortcomings: 

1. Insufficient sensitivity in low reflectivity hazard conditions, such as: 

a. Super-cooled, high liquid water content drizzle and clouds. 

b. High ice water content, high altitude clouds. 

2. Inability to correct for rain attenuation, resulting the under estimation of severe weather cores. 

3. Limited measurement range in highly attenuating heavy liquid precipitation. 

There are several possible design options that can be considered in the next generation of AWRs, but unfortunately 
there is no single design solution to improve detection in all the above listed hazard conditions. The following are 
the recommendations of this report, in order of increasing complexity and cost. 

Recommendation #1: Add Radiometer Measurement Capability to AWRs 

The addition of radiometric measurement capability to future AWRs should cost very little and can have potentially 
significant impact on detecting and avoiding inflight icing hazard, especially in super-cooled liquid clouds and 
drizzle that are often below the detection threshold of AWRs. The radiometric brightness temperature measurements 
can be used directly to identify the “coolest” direction with the least amount of liquid water path. Or if some 
assumptions can be made on drop size distribution (e.g., stratified precipitation), the radiometer measurements can 
be combined with the radar reflectivity profiles to quantify some cloud and precipitation parameters and thus the 
potential hazard. 



Recommendation #2: Support the Development of Dual-Polarized AWR Antennas to Enable Polarimetric 
Measurement Capability 

The use of polarimetry for characterizing precipitation (rain rate) and severe weather (hail detection) and correcting 
for attenuation has been well documented ([1-8]). The cost of upgrading the current AWRs to dual -polarized 
capability is relatively simple and inexpensive. The transmitted pulse has to be equally divided to the Vertical and 
Horizontally polarized antenna ports and a second receiver has to be added to the radar. The signal processing is 
computationally modest and polarimetric radar parameters, such as differential reflectivity (Zdr) and specific 
differential phase shift (Kdp) can be effectively used to correct for attenuation in moderate to heavy rain conditions. 
The only missing (enabling) technology is a dual-pol antenna, similar in quality to the current standard flat plate 
waveguide array. 


Recommendation #3: Develop Dual-Frequency (X/C-band) AWRs for Severe Weather Prone Regions 
In severe weather conditions X-band attenuation can be severe, resulting in underestimating dangerous weather 
cores, and worse, diminishing the received signal below the detectable threshold. By shifting to a lower frequency, 
attenuation quickly decreases and thus improves the probability of detection. A dual -frequency X/C-band system 
would have the higher sensitivity to smaller scatterers at X-band but the ability to penetrate highly attenuating (high 
rain rate) regions. 


1 Introduction 

Most AWRs (aircraft weather radars) operate at X-band (8-12 GHz) using magnetron tube or solid state amplifier 
transmitters and slotted waveguide array antennas. These X-band radars provide an optimal tradeoff between 
sensitivity, size, weight and cost, but are susceptible to severe attenuation in heavy rain and have limited or no 
ability to help pilots detect and quantify icing potential. This vulnerability of AWRs to attenuation and their inability 
to characterize and quantify precipitation has been a contributing factor in numerous aircraft incidents. 

Undetected peak convective cells, hidden by attenuation, led to the forced landing and resulting fatalities of 
Southern Airways flight 242 on April 4, 1977. If attenuation information was available from AWR, or if the radar 
signal was less attenuated, the peak convective region would have been correctly identified and likely avoided. 
Newer AWRs with Doppler capability attempt to work around this attenuation problem by measuring the radar 
signal spectrum width. Hazardous convective severe weather cores are almost always highly turbulent and thus 
cause a broadening of the Doppler signal spectrum. Furthermore spectrum width measurements are theoretically 
unbiased by attenuation, as long as the received signal is above the detection threshold. So the spectrum width 
technique can effectively detect convective severe weather, but attenuation can still diminish measurement range to 
impractically short distances. Since attenuation of atmospheric water droplets is approximately proportional to the 
square of the radar frequency, a potential solution to severe attenuation is multi-frequency AWR with a lower 
frequency channel. A second radar channel at C-band for example, would have a significantly lower rate of 
attenuation in heavy rain than at X-band, and consequently could better penetrate and more reliably detect severe 
weather. The multi-frequency radar technique is investigated in Section 4. Other methods to improve measurement 
sensitivity are discussed in Section 5. 

Some hazardous atmospheric conditions however do not have elevated turbulence or wind sheer, and radar 
backscatter intensity is not necessarily a good metric of hazard potential, such as freezing drizzle and rain. On 
October 31, 1993, undetected or under estimated freezing rain caused the crash of American Eagle flight 4184. 
Other icing related aircraft incidents are listed in http://www.skvbrarv.aero/index.php/In-Flight Icing . Icing hazard 
potential is related to the total volume of supercooled liquid along the flight path and drop size (larger drops can 
cause enhanced surface roughness and can have larger areas of impingement, sometimes beyond the ice protection 
system). A technique not yet utilized with AWRs for better measurement of severe weather precipitation or freezing 
rain is radar polarimetry. The quantitative measurement of rain with S-band (~ 3 GHz) polarimetric ground-based 
weather radars was successfully developed in the 1990s and by 2010 the polarimetric upgrade of the US network of 
159 ground-based weather radars began and is scheduled to be completed by the end of 2013 [1][2]. Moreover, new 
ground based radars are almost never built without polarimetric measurement capability. The initial low-frequency 
polarimetric rain retrieval techniques have been modified and validated for C and X-band [10] [15], but have not yet 
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been implemented with AWRs. One of the recommendations of this report is to upgrade future AWRs with 
polarimetric measurement capability for better characterization of rain in convective severe thunderstorms and in 
freezing rain. Weather radar polarimetry and its implementation with AWRs are discussed in Section 3. 

Supercooled clouds and drizzle can contain a large amount of liquid water comprised of small droplets, making 
them undetectable by most weather radars, yet still an icing threat to smaller aircraft. Higher frequency radars 
have better sensitivity to small droplets but also have a higher rate of attenuation in high liquid water 
conditions, making them actually less sensitive at the ranges that are needed from an aircraft, than the current 
X-band AWRs (see Section 0). Radiometers however have been effectively used to measure drizzle and cloud 
liquid for many decades and AWRs already have a receiver that can be used for brightness temperature 
measurements. This “passive”, receive -only measurement concept involves scanning ahead of the aircraft in 
azimuth while pointing in elevation above the horizon. The observed brightness temperature (or measured 
noise power) is the attenuation weighted, physical temperature along the antenna beam can be approximated as 
(see Eq. (5)) 

( 1 ) ^ \ 

where 

Cf j is the total atmospheric attenuation along the antenna beam in dB, and 

is the effective physical temperature of the attenuating cloud or precipitation layer. 

Since microwave attenuation due to small (relative to the wavelength) water droplets is directly proportional to 
the total liquid water, when the total liquid water along the beam is low, very low brightness temperature 
would be observed. If the total attenuation is high, the radiometer measurement will saturate, observing the 
physical temperature of the cloud or precipitation. But the total path attenuation can be controlled to some 
degree by pointing the beam up in elevation, to cut through a highly attenuating region at a steeper angle, 
thereby reducing the attenuating distance, and thus reduce the chance of saturation (Figure 1). If the 
propagation length through the cloud or precipitation can be approximated, then the average liquid density can 
be also estimated. But even without any knowledge of the propagation distance, by scanning the beam in 
azimuth at a constant elevation angle, the “coolest” direction can be found that corresponds to the safest path 
with the least amount of total liquid (Figure 2). This radiometric measurement technique is presented and 
analyzed next, in Section 0. 
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2 Radiometry 


The addition of radiometric measurement capability to future AWRs should cost very little and can have a 
potentially significant impact on detecting and avoiding inflight icing hazard, especially in super-cooled liquid 
clouds and drizzle that are often below the detection threshold of AWRs. The radiometric brightness 
temperature measurements can be used directly to identify the “coolest” direction with the least amount of 
liquid water path. Or if some assumptions can be made on drop size distribution, like in stratified precipitation, 
the radiometer measurements can be combined with the radar reflectivity profiles to quantify some cloud and 
precipitation parameters and thus the potential hazard. 


Observed Brightness 
Temperature 



Figure 1. Side view of an aircraft radiometer measuring a cloudy scene. As long as the antenna is 
pointed slightly above the horizon, the background is the cold space, making the observed 
brightness temperature highly sensitive to attenuation and therefore to the total liquid along the 

antenna beam. 
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Figure 2. Top view of an aircraft radiometer measuring the brightness temperature of a cloud-filled 
scene. The radiometer can scan in azimuth, slightly above the horizon, to look for the “coolest” 
direction with the least amount of total (supercooled) liquid. 

Radiometric Brightness Temperature Seen From Aircraft 

Figure 3 shows an aircraft flying at altitude. A, above the earth. Located in the nose of the aircraft is a 
radiometer that senses the brightness temperature, Tg in degree Kelvin (K), of Planck emission from the 
atmosphere and from objects in front of the aircraft. In the direction defined by vector, P, the radiometric 
brightness temperature, Tg, of the atmosphere can be computed as [3], 


where, 

6p = (see Figure 3), 

6g = 6p — 90° = elevation angle associated with P, 

P = distance from the aircraft along the vector, P, 

Pm = maximum distance along P thru the atmosphere (see text below), Kp(P) = Kp{P) + Ky(P) + Kq{P), 
total extinction (Np/km) at point P, 

Kp{P) = extinction (Np/km) due to liquid water droplets (e.g., clouds and rain), 

Kv(.P) = extinction (Np/km) due to water vapor, 

KoiP) = extinction (Np/km) due to oxygen, 

Tat(P) = air physical temperature (K), 

Pbg = background temperature (K) (see text below), 
and. 


( 2 ) 


TB = j^^Kr(P)-TAT{P)-e-^^^P^dP + T^Ge-r 
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( 3 ) 


T(a,fc) = j^Kriy) dy 
Tr = T{Q.PM) = dP 


(4) 



Figure 3. Geometry of aircraft flying at altitude, A, above the Earth. 


Note that Tj (Np) is the total atmospheric attenuation along the path defined hy P. Sometimes, (2) is 
approximated as [3], 


(5) Tgff[l—e + Tgg ■ e , 

where in degree K, is an effective physical temperature for the atmosphere. 

The integration in (2) along the path defined by P generally intersects many layers of the spherically-stratified 
Earth atmosphere. When P is tilted up from horizontal, this integration ends at the top of the atmosphere, and 
Pm is the distance from the aircraft (along P) to the top of the atmosphere. However, when P is tilted down 
from horizontal, the integration path can end either at the top of the atmosphere or on the Earth surface. When 
the integration path ends at the top of the atmosphere then is equal to the cosmic background temperature 
(~2.7 K). When the integration path ends on the Earth surface then is equal to the Earth brightness 
temperature, and P^ is the distance from the aircraft (along P) to the Earth surface. Earth brightness 
temperature varies over roughly 100-300 K, depending on surface composition, observation angle, and 
radiometer measurement frequency/polarization. 

Clearly, the proposed radiometer measurements should be conducted with the cold 2.7 K cosmos as the 
background to be sensitive to atmospheric path attenuation. So the antenna has to be pointed above the 
horizon. But, by how much? The next section (0) formulates a model for the gain pattern of a typical antenna 
and presents an example of measured earth surface bias as a function of elevation angle with a 1’ diameter 
AWR antenna in section 0. 
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The Radiometer Measurement 


The radiometer measurement is called antenna brightness temperature, (K). It is defined [3] as a spatially 
weighted average of Tg with weights derived from the antenna power-receive pattern, Here, the following 
calculation for will be used. 


(6) = /o "C (0 - f , </> - 7 t) ■ Ts da 

where, dfl = sin(0) cl6 dcp, dfl = 1, and (0, <p) are spherical coordinate system angles associated 

with the axes shown in Figure 3. Here, it is also assumed that. 


(7) 


Ga oc 


r sinCoii) 

[ £ti 


sinCttz) 1 ^ „ 


0-^ 


0.36 Og 


, «2 = 


(p—n 
0.36 Ob 
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where 6g (rad) is the antenna half-power beamwidth. 


Antenna beamwidth causes the radiometer field-of-view (FOV) to be a conical region. For the airborne 
radiometers shown in Figure 4, this conical FOV is depicted by the red -colored lines. Note that FOV for the 
radiometer on aircraft #1 makes visible the Tg of trees and soil on the Earth’s surface. For the radiometer on 
aircraft #2, the FOV does not contain the Earth but does make visible the sun, aircraft #3, and radar 
transmissions from aircraft #3. Here, the Tg associated with Earth, sun, and other aircraft are considered 
interference, since only atmospheric Tg is of interest. As suggested by Figure 4, there are strategies for 
avoiding this interference. For example, flying higher or pointing the radiometer antenna above horizontal can 
reduce interference from the Earth and from man-made transmissions (radars) located on the Earth. Solar 
interference can be avoided by pointing the radiometer antenna away from the sun. Software is available for 
determining sun location and can be used when the sun is obscured by clouds. Hence, with knowledge of 
aircraft altitude, heading, and pitch/roll, much of the radiometer interference can be greatly reduced via 
repointing of the radiometer antenna or mitigated by flagging the data. Finally, radiometers themselves can 
determine when the Ta measurement is corrupted by man-made signals. These interference detection 
techniques are described in [4] and [5]. 



Figure 4. Field-of-view (red boundaries) for radiometers aboard three aircraft. 
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Radiometer for Severe Weather Avoidance 


With the background of the previous sections in-hand, the radiometer’s role in severe weather avoidance can 
now be explained. First, note that radars can sometimes fail to detect severe weather because of intervening 
water-laden atmosphere that is both heavily attenuating (i.e., large Kj) and weakly reflective. This radar 
problem is shown pictorially in Figure 5 where the heavily attenuating atmosphere and severe weather (shown 
in red) has obscured radar returns from areas further in range (indicated by pink arrows). Flowever, as 
described below, the radiometer views the same scene much differently than the radar and would probably not 
be “blind” to the areas indicated by pink arrows. 
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Figure 5. Radar return showing two severe weather cells (red) that obscure returns from areas 

further in range (pink arrows). 



Note from (2) that when the atmosphere has no attenuation (i.e., Kj- ~ 0 in (2) or Tp ~ 0 in (5)) then Tb ~ Tbg- 
That is, without interference, the radiometer sees the cold background (T^g ~ 2.7 K) of outer space. Flowever, 
when the atmosphere has high attenuation (i.e., Kj- ~ oo in (2) or Tj- ~ oo in (4)), then Tg ~ T^ff ~ 250 K. As 
long as Tg lies between Tbg and Teff, then the radiometer measurement is said to be “not saturated” and a 
value of Tj can, in principle, be derived from Tg using (5). As shown in the sections to follow, an X-band 
radiometer measurement is seldom saturated in the proposed configuration. 

Radiometers cannot replace radars as a tool for severe weather detection, since radiometers cannot directly 
provide weather-distance (i.e., range) information. The following discussion elaborates on this point. Note 
from (4) that the value of the total attenuation (Tj-) derived from the radiometer measurement is the total or 
integrated value along the entire atmospheric path defined by P in Figure 3. Information in addition to the 
radiometer measurement is needed to determine whether the total attenuation (Tj-) results from a very deep 
weather cell having low attenuation rate (small Kj) or a thin weather cell having high Kj. Severe weather 
(e.g., heavy rain) is more strongly correlated with high Kj than high Tj. Flence, for severe weather detection, 
it is more useful to know how Kj varies along P rather than simply knowing Tj. The range information 
provided by radar measurements can be used to transform radiometer measurements of Tj into values of Kj. 
A novel method of deriving rough estimates of Kj using only the radiometer data is discussed in section 0. 
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However, should the radar fail, a novel all -radiometer severe weather avoidance strategy is possible. This 
strategy is implied hy Figure 2 and consists of scanning the area forward of the aircraft to find the region of 
minimum observed brightness temperature. The aircraft then flies towards this region since minimum 
brightness temperature implies minimum path attenuation as well as minimum liquid path (integrated rain rate) 
and thus (hopefully) the safest direction. 


Atmospheric Attenuation Details 

It appears most cost effective to upgrade existing X-band AWRs with either or both a radiometer mode and a 
dual-linear polarization (V and H) antenna. Hence, the focus here is on atmospheric attenuation at 9.3GHz for 
both V and H polarization. Before starting, we make clear the distinction between total atmospheric 
attenuation, in units of Nepers (Np) or decibel (dB), and atmospheric extinction (Np/km or dB/km). As 
shown in (3), the integral over distance of extinction produces attenuation. The units of extinction can be 
either Np/km or dB/km, where the conversion is 4.34 dB/km = 1.0 Np/km 

Plotted in Figure 6. is the calculated value of extinction caused by water vapor and oxygen in the US Standard 
Atmosphere (i.e., Ky + Kq in (2)) [3]. This calculation depends on altitude above the Earth because the density 
of both water vapor and oxygen thins with altitude and, to a lesser extent, because physical temperature at 
altitude has an effect. 

At the frequencies of interest in this report, dry-atmosphere extinction (i.e., Ky + Kq ) is typically much 
smaller than the extinction caused by rain and clouds. Rain and clouds can both be characterized as collections 
of liquid water droplets. The basis for extinction, caused by liquid water comes from Mie theory [3] [6], 
which is considered exact when the droplets have spherical shape. The spherical shape approximation is 
accurate for small droplets and the smallest of these remain suspended in the air as clouds. However, as shown 
in Figure 7, the larger droplets will fall through the air at significant velocity as rain. The air resistance at these 
fall velocities can change the droplet shape from spherical to oblate-spheroid. When the drop diameter nears 
0.005 m, the fall velocity can cause sufficient air resistance/turbulence to break the droplet apart [7]. The 
oblate-spheroid shape of the droplets causes the extinction of rain to be polarization dependent. Convenient 
formulas that take this polarization dependence into account are given by [8]. These formulas were used to 
calculate the rain extinction at both V-pol and H-pol experienced by a 9.3 GHz wave traveling in the horizontal 
direction thru rain. Results of these calculations are plotted in Figure 6.. 
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Extinction (dB/km) 

Figure 6. Extinction at 9.3 GHz for rain and for the rain-free US Standard Atmosphere 



Figure 7. Raindrop shape depends on size and fall velocity. Drops larger than 5mm 
in diameter often break into smaller drops [7]. 

Simulated Radiometer Measurements in Severe Weather 


Software was developed to compute values of that are seen by an airborne radiometer. The software uses 
formulas from sections 0 and 0 and can account for the following factors: 

• radiometer operating frequency 

• antenna beamwidth and gain pattern 

• antenna pointing (elevation angle) 

• aircraft flight altitude 

• round Earth effects 
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• Earth surface brightness temperature 

• atmospheric physical temperature, pressure, water vapor (all functions of altitude) 

• rain rate 

• rain height 

• rain storm radius 

• aircraft location within the rain storm 

Here, calculations from the software are used to determine what conditions (if any) cause the radiometer 
measurement to saturate at as discussed in section 0. The following notes apply to these calculations. 

1. The radiometer operating frequency is 9.3 GHz, which is the operating frequency of many airborne 
weather radar (AWR). 

2. Antenna beamwidth of Og = 7.5° = 0.13 rad is used and is typical for a medium sized AWR antenna. 

3. The rain height is 3 — 4 km, which is typical for many mid-latitude storms. Rain height is also 
referred to as the height of the melting/freezing layer. Above this height, the water is in its frozen 
state. Frozen water causes much less attenuation than the liquid water located below the rain height 

[3]. 

Results of our first calculations are plotted in Figure 8 and are for a small severe weather cell, similar in size to 
a small thunderstorm. In particular, this weather cell has a rain height of 3 km and the rain-filled region of the 
cell extends to a range of 12 km in front of the radiometer/aircraft. The Figure 8 (a) data results when the 
radiometer antenna is pointed straight ahead (i.e., 6^ = 0°). The three curves in this plot correspond to three 
separate flight altitudes, with the lower (1 km) and upper (5 km) altitudes being below and above the rain 
height respectively. The Figure 8 (b) data results when the antenna is pointing up by about one-half of the 
antenna beamwidth (i.e., 6^ = 4°). For these calculations, Teff ~ 280 K and represents the saturation 
value. Note that Tj^ saturation is observed when the radiometer is flying below the rain height and only when 
rain rates are extreme (> 60 mm/hr). Also note the differences between Figure 8 (a) and Figure 8 (b) at low 
rain rates, which are due to differences in 6^. That is, for the slight up-tilt of the antenna in Figure 8 (b), Tj^ 
becomes less biased (corrupted) by the Earth brightness temperature (here assumed to be 288 K). 

Results of our second calculations are plotted in Figure 9 and are for a large thunderstorm. For this storm, rain 
height is 4 km and the rain-filled region of the cell extends to a range of 35 km in front of the 
radiometer/aircraft. As in Figure 8, Tgff ~ 280 K and Earth brightness temperature is 288 K. Unlike Figure 
8, the flight altitudes of 1 km and 3 km are now both below the rain height. Also note that, due to the larger 
extent of the storm, Tj^ saturation occurs when rain rates are very-heavy (> 30 mm/hr) instead of extreme. 
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Antenna Elevation = +4 deg 



0.1 1.0 10.0 100.0 1000.0 
Rain Rate (mm/hr) 

(b) 


Figure 8. Computed brightness temperature of a SMALL thunderstorm for antenna elevation 
angles of (a) Odeg and (b) 4deg and for flight altitudes of 1km, 3km, and 5km. The rain height is 3 

km. 


12 




Antenna Elevation = 0 deg 



0.1 1.0 10.0 100.0 1000.0 
Rain Rate (mm/hr) 

(a) 


Antenna Elevation = +4 deg 



0.1 1.0 10.0 100.0 1000.0 
Rain Rate (mm/hr) 

(b) 

Figure 9. Computed brightness temperature of a LARGE thunderstorm for antenna elevation 
angles of (a) Odeg and (b) 4deg and for flight altitudes of 1km, 3km, and 5km. The rain height is 4 

km. 

Retrieving Total Atmospheric Attenuation from Radiometer Data 

In (5), the substitution, Lj = is used to give (9) where Lj is interpreted as total loss of the atmosphere. 


( 8 ) Ts ^ Teff ■[!- It] + Tbc It 

Straightforward manipulation of (8) gives the following solution for Lj. 


( 9 ) 


^ _ Tg-T^ff 

^ TsG-Teff 
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The precision, AL, associated with a retrieval of Lj from Tg using (9) depends on the measurement precision 
(NEDT), AT, for Tg. See Section 0, for further discussion on AT. The relationship between AL and AT can he 
derived from (8) as. 


( 10 ) 


AL = 


AT 
\ 9Tb \ 
\ dL-p \ 


AT 

Teff-Tsc 


Using Teff = 280K, Tbq = 2.7K, and two values of AT (NEDT), AL can was calculated from (10) and plotted 
in Eigure 10. 



Figure 10. Atmospheric attenuation retrieval precision caused by radiometer NEDT. 


However, AL can also he affected by errors in Teff- Eor example, using a value of Teff that is either too high 
or too low by 1% (~3 K) will bias the retrieved value of Lj as shown in Eigure 11. When reviewing Eigure 
11, note that a 100% error in Lj- is also a 3 dB error in Lj. Thus, for example, when actual Lj- w 20 dB, a 
100% bias in Lj implies either 23 dB or 17 dB. In principle, when the radiometer antenna is pointing at small 
elevation angles, the physical temperature of air near the radiometer/aircraft should serve as an accurate 
estimate for Teff- 
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Figure 11. Atmospheric attenuation retrieval bias caused by error in Teff- 


Retrieving Rain-Rate and Rain-Depth from Radiometer Data 

As discussed in section 0, radiometer measurements are typically used to estimate total path attenuation, 
Tj = — In(LT-). However, estimates of attenuation rate, Kj are more useful for severe weather detection. 
Knowing Kj is equivalent to knowing both rain rate, R (mm/hr), and rain depth, (km). Here, a novel 
method is presented for estimating R and from only measurements of the radiometer. 

Using the relationship in (5), the V-pol and H-pol measurements of a radiometer that is viewing rain of depth, 
du, can he approximately calculated as, 

(11) Tb(p) « . 

Where p is the polarization, commonly V (vertical) or H (horizontal). 

From [8], one finds that at 9.3GHz, 


(12) K,(P)= Kp- R‘^p , p = (V,H), 

where, V — pol: (KY,ay) = (0.00792,1.265) 
and H — pol: (k^, a^) = (0.00881, 1.296). 

Hence, Tg (V) and fg (H) represent two equations of the form shown by (1 1) that contain two unknowns, R and 
djj. The closed-form solution for R is given by. 


(13) 


R = 



ln\ 

'^eff Y 

ay-an 

KH 



Kv 

2n| 

'^eff \ 



TefrTBm)\ 
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Once a value of R is computed from (13), it can be substituted into (11) to solve for d^. Note that the 
condition, ay ^ a^, makes the solution for both R and dg possible. When ay = = aQ, then there is 

effectively only the one unknown, dg ■ and separate solutions for dg and R are not possible. 


Further analysis shows that the solutions for d^ and R from 9.3 GHz radiometer data described above are quite 
sensitive to AT. Retrieval error in dg and R for AT = 1 K in both Tg(V) and ^^(H) was computed and is 
shown in Figure 12. As seen in these plots, retrieval error depends on the actual values for dg and R. 




Figure 12. Retrieval error for (a) rain rate and (b) rain-depth when raw radiometer measurement 

error is 1.0 K. 
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Brassboard Radiometer Description 

For experimental purposes, a simple radiometer (i.e., brassboard) was constructed. The radiometer is described by 
the block diagram in Figure 13 and by the photo in Figure 14. The radiometer antenna was made by Cobham (part# 
460386) and is commonly used by aircraft weather radar. Receiver components were housed in a Styrofoam 
enclosure to help stabilize receiver physical temperature and reduce radiometer calibration drift. Data acquisition 
from the radiometer utilized an Agilent Technologies spectrum analyzer, part# E4440A. The radiometer was 
calibrated using commercially-available calibrated noise sources and microwave absorbing material. Calibration 
data is shown in Figure 15. 

Error! Reference source not found.Brassboard radiometer receiver. 



Figure 13. Block diagram of the brassboard radiometer. 
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Figure 14. Brassboard radiometer. 
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Figure 15. Calibration curve for the brassboard radiometer. 


Brassboard Radiometer Experiment 

The brassboard radiometer described in the previous section was used to validate our theoretical description of 
the radiometer measurement (Tj^) given in sections 0 and 0. More specifically, data was collected from the 
ground-based brassboard radiometer as it was tilted to scan from Earth to the sky across the horizon. The 
effective horizon during this experiment was the top of the wide building shown in Figure 14. Measured and 
calculated values of for this experiment are shown in Figure 16. First, note that the horizon for this 
experiment was measured to be at an elevation angle of 5°. Thus, when the radiometer antenna was scanned 
to 6£- = 5° then half of its antenna beam observes the Earth and half observes clear sky. Hence, the measured 
Tyi at 6£- = 5° is seen to be about half-way between the clear sky value of ~3 K and the Earth value of 

18 


~250 K. Calculations plotted as a solid line in Figure 16 used the actual radiometer antenna beamwidth of 
= 7.5° = 0.13 rad. These calculated values of agree well with the measured data, giving us confidence 
with our modeling. The dashed-line calculations shown in Figure 16 used 6^ = 0° and can be interpreted as 
the Tg -profile of the scene. As discussed in Section 0, the brassboard radiometer Tj^ measurements represent 
an antenna-beam weighted average of this Tg -profile. These results indicate that an AWR antenna only needs 
to be pointed up in elevation by less than the beamwidth to prevent loss of sensitivity from earth surface 
radiation. 



Figure 16. Measured brightness temperature as a function of elevation angle. This shows that a 7.5 
deg beamwidth (1’ diameter) antenna has to be pointed in elevation a minimum of about 5 deg 
above the horizon to prevent significant earth surface radiation interference. A larger antenna with 
narrower beamwidth could be pointed even closer to the horizon. 


Combined Radar/Radiometer Retrieval 


In some cases co-located radar and radiometer measurements can be made. Flere “co-located” emphasizes that 
a similar range interval is sampled and detected by the radar and that the same range interval contributes most 
of the path attenuation and thus the radiometer measured brightness temperature. The next analysis shows how 
such co-located radar and radiometer measurements might be combined to better characterize liquid 
precipitation. 

The general radiometric transfer equation for the radiometer observed brightness temperature is given by [24], 
Vol. 3, Ch. 17, pp. 1331, Eq. 17.80 as 


(14) 


CO 

T,=\K,(p)r,,(p)e 


-\Kj{P’)dP' 


dP + 


00 

-\Kr{P)dP 
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where 


^BG 


is the background radiation , 


f) f) 

^ is the Zenith angle ( ^ = 90 — > pointing towards the horizon), 


K 


^ is the total atmospheric absorption coefficient in units of nepers per unit distance at distance P from the 
radar, and 


is the temperature of the atmosphere at P. 

In a stratified cloud/precipitation, where uniform physical temperature and drop size distribution can be 
assumed, the observed brightness temperature can be approximated with the following much simpler 
expression [24] Vol. 3, Ch. 17, pp. 1331, Eq. 17.83 : 


( 15 ) 

where 


T =T e 


-k„AP 


7 ’ 1 


^47 is the effective (or mean) radiating atmospheric temperature, and 

is the extent of the cloud/precipitation layer along the antenna beam direction . 




The exponential attenuation factor, ^ through the cloud-precipitation along the antenna beam can be 

/ = ^ ° 

defined as the total loss factor, ‘ , and since the comic background radiation is a negligible 2.7 K, 

the observed brightness temperature can be simply approximated as 


( 16 ) 




So if the mean cloud/precipitation physical temperature can be estimated from in situ measurements and 
extrapolated into the cloud/precipitation layer, then from the observed brightness temperature, the total loss 
factor can be calculated. Lj can be also expressed as a function of the more commonly used attenuation 
coefficient, a in units of dB/km, where Ka (nepers/km) = 0.23 « (db/km) 

^-0.23aAP 

f 

( 17 ) 


Lj, =e 


From the radar measurements, the cloud/precipitation path length can be directly measured and then the 
path attenuation rate a can be calculated. 

The attenuation rate is directly proportional to the liquid water content for hydrometeors that are much smaller 


n27tr/ 


a 


<0.5 


than the radar wavelength, Rayleigh scatterers (when 
cloud/precipitation drop radius and X is the radar wavelength) [3] : 

a = m ,/‘ exp [- 6.866(l + 4.5 x lO'^r)] 


, where n is the index of refraction, r is the 


( 18 ) 

where 


Tny = liquid water content (g/m^). 
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/ is the radar frequency and T is the cloud/precipitation temperature. 


So from the radiometer observed brightness temperature and radar measured cloud/precipitation path length, 
the attenuation rate and then mean liquid water content can be estimated. 


m„ 


(19) 



— ND^ 

6 


where 


Di is the i* cloud/precipitation scatterer diameter, Dy is an equivalent volume diameter and Ny, the equivalent 
volume number density. 


The radar received power for Rayleigh scatterers is proportional to the sixth moment of the drop size 
distribution (see section 4.1.2, equation (68 ): 


( 20 ) 


^ PfAn^cT ^ 

' 32r" h 


n' 




r 


And similarly, defining an equivalent volume diameter, Dy, and number density, Ny, the radar measured power 
can be written simply as 


Pr = 


P^An\cT7t K 


( 21 ) 


2 o4 


32rU' 


■NyO; 


By combining the radiometer measured liquid water content with the radar received backscattered power, the 
equivalent volume diameter and number density can be estimated: 


( 22 ) 


3P^AnlcT7i^ K 
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SO 

( 23 ) 

and 

(24) 
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A/ - 3 I 


P 16 rU 


my3P^An^cT7T K 
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A measure of drop size and liquid water content are important parameters for characterizing the hazard 
potential of super-cooled precipitation. Higher liquid water content clearly presents an increased hazard 
potential, because there is more liquid available to freeze on the aircraft wing. But drop size is also important, 
because larger drops shatter on impact and can cause enhanced surface roughness and have larger areas of 
impingement, sometimes beyond the ice protection system. 
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Adding Radiometer Mode to a Weather Radar System 

New -model, digital-output, AWR (e.g. Rockwell RTA-4100) can cost $150,000. Old-model, digital -output 
AWR (e.g., Honeywell PRIMUS -440) cost $80,000 for a new system and used systems cost about $30,000. 
Digital-output AWR have mostly replaced the older analog-output AWR. However, some analog-output 
AWR are still utilized and have a used price of roughly $5,000-$15,000. Major manufacturers of aircraft 
weather radar (AWR) are currently: 

• Rockwell-Collins 

• Honeywell 

• Bendix-King 

• Telephonies 

• Garmin 

To permit planning for the addition of a radiometer mode to an AWR, technical information was requested 
from some of these AWR manufacturers. However, the contacted manufacturers required a non-disclosure 
agreement from us before providing any technical information. This disclosure agreement seemed 
inappropriate at this early planning stage. So, detailed AWR technical information was not obtained. 
However, as shown below, many aspects of joint AWR/radiometer operation can still be investigated. 

The method of using a radiometer to assist the AWR with severe weather avoidance was described in sections 
0 and 0. For this application, moderate measurement accuracy (~5 K) is required from the radiometer. In 
principle, existing radar circuitry could be used to make the radiometer measurements and achieve this 
accuracy after some simple hardware modifications. These hardware modifications consist of adding a switch 
and calibrated noise source to the AWR front-end electronics as shown in Figure 17. Note that noise source is 
used as a warm load calibration unpowered, and as a hot load when the power to the noise source is turned ON. 
The radar receiver is then reprogrammed to occasionally toggle both the new switch, S2, and the noise source 
power so that data can be collected from the noise source in both ON and OFF state. This noise source data is 
then used for radiometer calibration in the manner described by Figure 15. As shown in Figure 18, operation 
of the AWR is also modified so that the transmit (TX) pulse is absent or blanked immediately before 
radiometer data collection. Following the blanked-pulse, all the range -gate data samples collected over the 
pulse-repetition-interval (PRI) can be averaged and used as radiometer data. 
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Figure 17. Aircraft weather radar (a) before and (b) after radiometer mode modification. 
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Figure 18. Blanking the transmit (TX) pulses for the radar permits radiometer data collection. 


With this overview in hand, details of radiometer operation and performance can now be discussed. Precision, 
AT (K), of the radiometer measured 7^ can be estimated from [3], 


(25) 


_ Ta+ Trec 
y/BlF' tA 


where Trec (K) is the receiver noise temperature, Bip (Hz) is the bandwidth of the receiver, and (s) is the 
time interval over which the data is averaged. As shown in Table 1, modern AWR typically have K 

and Bif~S20,000 Hz. When Tj^ is near 280 K, calculation (25) shows that achieving AT~5 K will require 
tyi~0.025 s. Further, if the radar scans a 90° (azimuth) X 30°(elevation) sector forward of the aircraft with a 
7.5° wide antenna beam, then about (90/7.5) X (30/7.5) = 48 separate radiometer measurements are needed. 
Since each measurement takes 0.025 s, a total of 48 X 0.025 = 1.2 s will be required to perform a complete 
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scan. Thus, the aircraft pilot would have severe weather updates from the radiometer at a rate of about once a 
second. 


However, surely the radar receiver bandwidth could be opened up to at least 5 MHz, and modern data 
acquisition system can easily process data at that sample rate. So repeating the above calculation, the 
approximate update rate with this more realistic receiver bandwidth would decrease to 0.4 sec, yielding two 
scan updates per second. 


Table 1 Specifications of a modern aircraft weather radar. 


Design Parameters 

Derived Parameters 

Frequency (MHz) 

9500 

Wavelength (cm) 

3.16 

Pulse Width (ps) 

2 

Receiver Bandwidth (MHz) 

0.52 

Pulse Repeat Freq (kHz) 

6 

Maximum Range (km) 

25 

Rcvr Noise Figure (dB) 

3 

Rcvr Noise Temperature (K) 

290 


Radiometer Summary and Conclusions 

Adding a radiometer mode to conventional aircraft weather radar (AWR) will improve hazardous weather 
detection/avoidance by overcoming some radar deficiencies (e.g., insufficient sensitivity in liquid cloud and 
drizzle conditions, severe weather blind-spots) and assisting the radar to determine atmospheric attenuation. 

Radiometer measurements can be processed into estimates of total atmospheric attenuation (dB). However, as 
discussed in section 0, atmospheric extinction (dB/km) is a more useful indicator of severe weather than 
atmospheric total attenuation. Radiometer measurements can be combined with radar measurements to derive 
atmospheric extinction. A novel method for roughly estimating extinction using only radiometer data was 
presented in Section 0. However, it is generally more accurate to use the combined measurements from radar 
and radiometer to estimate extinction. 

A novel radiometer-only strategy for severe weather avoidance was presented in section 0. This all-radiometer 
strategy consists of flying towards the area of least atmospheric attenuation and should help overcome radar 
blind-spots. It was calculated that radiometer blind-spots (i.e., saturation) occur only under the unusual 
conditions of extreme rain rates that extend more than 12 km in front of the aircraft. Tilting the antenna further 
up in elevation reduces the chance of saturation even in severe conditions.. 

Radiometer measurements can be contaminated by interference from the sun. Earth, other aircraft, and man- 
made transmissions (e.g., radar). However, most of this interference can be avoided or flagged by monitoring 
and controlling the antenna pointing. Much of the remaining interference can be self-detected by the 
radiometer. 

A brassboard radiometer was constructed and used to collect data during a horizon-scanning experiment. 
Radiometeric measurements and calculations from this experiment show excellent agreement and give 
confidence in our modeling. Furthemore, the measurements demonstrated that an AWR antenna only has to be 
pointed above the horizon by less than the antenna beamwidth to prevent loss of sensitivity due to earth surface 
radiation. 

Although we have not been able to obtain detailed technical information from vendors of AWR, planning the 
addition of a radiometer mode to a typical X-band AWR was investigated. Calculations show that severe 
weather updates from this radiometer can be obtained at a minimum rate of about 2 Hz for 48 separate cells in 
90° (azimuth) x 30°(elevation) sector forward of the aircraft. 
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3 Polarimetry 


Review of Weather Radar Polarimetry 

Current aircraft weather radars collect radar reflectivity (power) and Doppler velocity mean and standard 
deviation (changes in signal phase) at horizontal linear polarization. Hazardous weather conditions are 
identified hy high reflectivity and/or hy elevated velocity standard deviation (high turbulence or sheer). 
Reflectivity however is unreliable in severe weather conditions, because rain attenuation can mask severe 
weather cores and icing hazard (freezing rain) for example is not necessarily accompanied by elevated 
spectrum width. Dual-polarization radar can measure additional radar parameters that can be used to 
correct for signal attenuation in rain and can identify mixed phase precipitation such as hail (and the 
melting layer). 

The polarization of a radar system refers to the orientation of the electric field oscillations, perpendicular 
to the transverse wave’s direction of propagation. Current weather radars operate using horizontally 
polarized antennas, so when the radar is pointed towards the horizon, the electric field oscillates in the 
horizontal plane as shown in the first (left) graphics of Figure 19. The transmitted field of dual- 
polarization radar, shown in the second (right) graphics of Figure 19, uses both horizontal and vertical 
polarized fields. The horizontally polarized filed scattering most influenced by the horizontal cross 
section of particles, while the vertical polarization “sees” the vertical dimension. Small raindrops are 
close to perfect spheres due to the surface tension, but larger deform and often modeled as oblate 
spheroid. This deformation is illustrated in Figure 20. When only small hydrometeors are present in a 
range cell of polarimetric radar, the backscattering, attenuation and phase propagation is approximately 
same at H and V polarization. But when larger raindrops are present, the H-pol radar pulse will reflect 
stronger, attenuate more and phase propagation will be slower than V-pol. Consequently, the received 
power ratio (Z^r) and the relative phase {<i>pd) between H and V-polarized received signal contain 
information on hydrometeor size. The phase center of mixed phase and non-water particles also different 
at V and H-polarization, so mixed phase precipitation (hail and melting band) and ground debris can be 
effectively detected and differentiated from rain and drizzle using H and V signal correlation (phv)- The 
following sections describe the polarimetric radar parameters commonly used for characterizing 
precipitation. 




Figure 19. Dual polarization technology gives additional information on hydrometeor size, shape 
and composition to allow the estimation of rain rate, signal attenuation and the detection of hail. 

Source: NOAA 
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Figure 20. Raindrop shape primarily depends on size. Larger raindrops also have significant shape 

oscillations and often break apart. 


3.1.1 DIFFERENTIAL PHASE SHIFT (0dp) 

Odp is the difference in the returned phase of the vertical and horizontal returns. The difference in return 
phase is greater from areas of large ohlate spheroid drops because the horizontal pulses are slowed 
slightly more than the vertical pulses. This slight difference is enhanced hy range from the radar. It is 
generally regarded that Odp is best used in its range derivative form: Kdp. 


3.1.2 SPECIFIC DIFFERENTIAL PHASE SHIFT (KDP) 

Kdp is the range derivative of Odp in units of deg/km. It is generally immune to large hail and thus is 
considered a useful indicator of actual rain rates and thus attenuation. Larger values indicate higher rain 
rates and higher attenuation. Light rain containing smaller, non-oblate spheroids results in very little 
return phase differences. 
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Figure 21. Dependence of Rain Rate - Kdp relations on radar frequency, (from Park, S.-G., Maki, 
M., Iwanami, K., Bringi, V.N., Chandrasekar, V., 2005. Correction of radar reflectivity and 
differential reflectivity for rain attenuation at X band. Part I: theoretical and empirical basis. J. 

Atmos. Ocean. Tech. . 22, 1621-1632.). 


3.1.3 DIFFERENTIAL REFLECTIVITY (Zor) 

Zdr is ten times the log of the ratio of the horizontal and vertical reflectivity return from a sample 
volume. Large drops are ohlate spheroids and have increasingly larger returns from the horizontal pulses 
compared with the vertical pulses. For Zdr over 6 dB, the drop diameters are exceeding 8 mm. Drops with 
diameters under about 0.5 mm have Zdr about zero. Unfortunately, hail also has a Zdr of around zero. 
Generally, snow has a low Zdr independent of shape due to its small dielectric constant. Zdr is sensitive 
to drop size, but inversely to propagation through rain with larger raindrops. So the observed Zdr first has 
to be corrected for differential attenuation and then usually combined with Kjp (specific differential phase 
shift and/or Z (reflectivity) to estimate rain rate/attenuation. The linear reflectivity ratio is usually 
differentiated from the LOG ratio by using lower case “dr”: 

„„ =101ogio(z*) 


3.1.4 CORRELA TION COEFFICIENT (pav) 

Phv is a measure of the correlation of the returned power from the horizontal and vertical pulses. Typical 
meteorological values are near unity with a very small variation. In regions of non-uniform scatterers, 
such as mixed phase areas (i.e. - melting layers), the pnv is lower. Ground clutter and airborne debris also 
have lower values as shown in middle image in the top row of Figure 22. pev is rather used for scatterer 
identification than quantitative retrieval. 
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Figure 22. Example of X-band polarimetric radar data collected at low elevation angle in a 
tornado: Reflectivity (top-left), phv debris signature (top-middle), spectrum width (top-right), 
Doppler velocity folded (bottom-left and ) and un-folded (bottom-middle), Zdr (bottom-right). 
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Rain Rate Estimation from X-band Polarimetric Radar Parameters 


Rain rate is highly correlated with attenuation, so existing models for estimating rain rate from X-hand 
polarimetric radar parameters can he used to correct for attenuation. 



Figure 23. examples of k-R relationships established using Mie theory (spherical raindrops) for the 
K- (0.86 cm), Ka- (1.15 cm), X- (3.2 cm), C- (5.6 cm), and S-band (10 cm) wavelengths and the for 
raindrop temperatures T=0°C (dotted line), T=10°C (continuous line), and T=20°C (dashed line). 

From [12]. 


3.1.5 Rain Rate - Z Estimator 


(27) 

from [13] , 


R = 100z;f 


where R is rain rate in mm/hr, and 


Zeh is equivalent reflectivity factor at horizontal polarization, in units of mm®/m^. 
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3.1.6 

(28) 


Rain Rate - Kjp Estimator 


R = aK 


DP 


where 

Kjp is the specific differential phase shift in units of deg/km, and 
a = -16.5 and b-~ 0.75 hut slightly dependent on drop size distribution. 



R=17.0 K (Knp>0.1 deg km'^) (HMT-06) 
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Figure 24. R - Kdp power-law fit relations using size distributions collected in various field 

campaigns. From [14]. 

3.1. 7 Rain Rate - Z, Zdr Estimator 


R = 0.0039Z:^’Zjf ’ 

(from [18]) 
where 

Zjr is differential reflectivity in linear units; Zdr - lOlogio(Zrfr). 


30 


3.1.8 Rain Rate - Z, K^p Estimator 

R = a ■ 

(30) ‘‘P 

(from [19]) 

where the a, h, c, d coefficients are given in Table 2. 


Table 2. Coefficients of Equation (30) rain rate estimator (Kdp, Z, Zdr). 


Rain type 

a 

B 

C 

D 

Stratiform 

0.009 

-0.173 

0.103 

-0.653 

Convective 

0.198 

0.4405 

0.035 

-0.036 


3.1.9 Comparison of various rain rate estimators 

The various rain rate estimators are compared in Figure 25 for sensitivity to variability in drop size 
distribution. As expected, the reflectivity only estimator is very poor in conditions anything other than the 
lowest rain rates. The estimation improves, as more parameters are included with reflectivity. But only 
error free parameters are considered in these simulated results, while Z and Z± are susceptible to error due 
to unaccounted for attenuation and differential attenuation. Kdp on the other hand is a much easier to 
measure accurately, so the upper right rain rate estimator plot is closer to what can be realistically 
expected than the two lower row plots of the more complex models. 
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Figure 25. Scatter plots of rain rain calculated from measured drop size distribution and estimated 
from four types of rain rate estimators a) R(Z), b) R(Z, Zdr), c) R(Kdp), d) R(Kdp, Z, Zdr). From [19]. 

Dual Polarized Aircraft Weather Radar Antenna. 

The cost of upgrading current AWRs to dual-polarized capability is relatively inexpensive. The 
transmitted pulse has to be equally divided to the Vertical and Horizontally polarized antenna ports and a 
second receiver has to be added to the radar. The signal processing is well-documented and straight 
forwards and polarimetric radar parameters, such as differential reflectivity (Zjr) and specific differential 
phase shift (Kdp) can be effectively used to correct for attenuation in moderate to heavy rain conditions. 
The only missing {enabling) technology is a dual-pol antenna, similar in quality to the current standard 
flat plate waveguide array. A dual-pol slotted waveguide array antenna has been patented by Rockwell 
Collins: http://www. google.com/patents/US8Q98 1 89 but performance data was not available to assess 
potential polarimetric measurement accuracy, at the time of this report. 

Past aircraft weather radars used magnetron transmit sources that required waveguide antenna array 
technology to be able to transmit high peak power. Modem weather radars use lower peak power solid 
state transmitters, which opens up the possibility of using microstrip antenna arrays. Some attempts have 
been made to develop lightweight dual-pol X-band microstrip antennas. In 2001 a design study was 
conducted by Antenna Design Associates to develop a Multi-Frequency (X, Ka and W-band) and dual-pol 
at X-band aircraft weather radar antenna. It concluded that a two-layer aperture coupled microstrip array 
technology (shown in Figure 26) could be used to develop a 1’ to 2’ diameter, medium power (<10 kW 
peak power) flat, lightweight X-band aircraft weather radar antenna. Larger antennas with narrower 
beam-width would be also possible, but without much improvement in gain due to increasing losses, with 
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increasing antenna size, in the microstrip substrate. More recently, a small experimental dual-pol 
microstrip array (~5” diameter, 25 deg beam width, 19 dB gain) was developed by [26]. 



Figure 26. Dual-polarized X-band microstrip antenna array with interleaved coplanar feed (coax 
feed points are shown as solid circles; offset = 6.96 cm = 0.202D). This design also incorporates a 
0.2D V-H aperture offset for spaced antenna transverse wind measurement. A waveguide power 
divided network is also required in the rear of the antenna. 

3.1.10 Antenna requirements 

Current airborne weather radars use flat plate (or box -horn) array antennas, like shown in the left pane of 
Figure 27. These antennas are ideal for airborne radar applications due to their light weight, low profile 
and excellent electrical properties (low loss and low sidelobes), but are not available in dual -polarization 
version (no information is available on the Rockwell Collins patent of a dual-pol slotted waveguide array, 
http://www.google.com/patents/US8Q98189) . Microstrip antennas are also relatively light weight and 
have low profile, but can be developed with dual polarization capability, so they might be a potential 
candidate for future dual-polarized aircraft weather radars. The following discussion is on the required 
specifications of a dual-pol AWR antenna for hydrometeor classification and implementation of 
polarimetric rain retrieval algorithms. 


33 



Figure 27. Various size flat plate (box-horn) waveguide array antennas. 

An ideal weather radar antenna has low loss, low sidelohes, close to the ratio of wavelength ( ^ ) and 


antenna diameter (D), (rad) half power heamwidth (high gain for a given diameter) and if 

polarimetric, then the orthogonally polarized ports have high isolation, hut co-located beams. The high 
gain is important for sensitivity. The low sidelohes (>25 dB one way) are important to ensure confidence 
in the location of the source of the measurements. Sidelohe return from a strong reflectivity core can 
severely bias an observation. Fortunately most well designed antennas have at least 20-25 dB first 
sidelohe levels and sidelohes further from the main beam quickly drop. But more importantly, with single 
antenna radar, the antenna is used for transmission and reception, so interference from a precipitation core 
is attenuated by twice the sidelohe suppression level in dB. Consequently an antenna with a 25 dB 
sidelohe will suppress return from the direction of that sidelohe by 50 dB. Since most airborne weather 
radars rarely measure above 50 dB signal to noise ratio, a 25 dB first sidelohe level is usually sufficient. 

Co-located V and H-pol beams are essential for useful measurements, with correlations varying 
from 0.99 in light rain and drizzle to 0.90 in hail and as low as 0.80 in mixed phased precipitation 
(melting layer / bright band). Consequently the vertical and horizontally polarized beam patterns of a 
dual-pol antenna should be 99% co-located (overlapped). 

Quantifying the effect of polarization isolation on the measurement accuracy of polarimetric radar 
parameters and consequently on the quality of the retrieved rain rate or attenuation is more complicated. 
Wand and Chandrasekar investigated this problem [21], but only plots of the results were published, not 
the complete expressions. Consequently the effect of finite antenna isolation on rain and attenuation 
estimation required us to re-examine the problem: 

The received field of simultaneous V and H radar transmission can be described as 


( 31 ) 
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where 


Q 

^ is transmitter phase difference between the V and H channels, and [A] is the antenna distortion matrix, 

Q 

accounting for antenna polarization isolation, /, V-channel phase shift and coupling phase shifts and 


VH . 


w= 


(32) 


1 
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j^vi 


j^H 


,j^\n 


[S] is the target scattering matrix: 


[ 5 ] = 


(33) 


S' S' 

S' S' 


'F 

and [P] is the propagation distortion matrix, accounting for differential phase shift ( 


w= 


(34) 


1 0 


0 e ' 


The resulting received voltages from the H and V antenna ports are: 


(35) 


R. = 


SvH^ 




Sfjfj (l + a/7 )+ Sjjye ^ )+ 

41 (l + 41e^^"''e ^^^ )+ Syye~^'^“'’ 41e^^^'' + 41e^^''^ ) 


and 
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These expressions can he simplified hy eliminating the second order terms with 

s^41 , 


{Jlf Sy„jl 


and 


/ \ ^ / 

(37) Rh = Shh (l + a/7 4^^ )+ Sfjye^^ 4^^" 4^^ + Syye~^'^“’’ 41 4^^^ 


(38) Ry = SyyC 4^^ + 41 4^™ )+ Sy^^e ^ 41 4^^'' 4^^ 
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The resulting becomes a simplified expression for the observed radar reflectivity factor 2 by assuming 
that co-pol scattering is uncorrelated with cross-pol backscatter: 


/ \ ^ / 

(39) (l + a/7 )+ + Syye~^'^‘"’ 4l 

/ \ ^ / 

(40) Ry = Syye~^'^‘“’e^^'^ + a/7 )+ Sy^e~^ + S^h V7 


(41) ^ 


1 + JJ{e ^ + e ^ + 


Pm. 

V V ^dr 


J_ ^jOvv ^jSt ^j^VH _|_ ^+j'^dp ^-iOvv ^-j0T ^-j^VH 


) 


and the complex exponentials combined to simplify: 




(42) 


1 + 2a/7 COs(6j,y +0J.)+ 2p^y — COs(T^^ -6yy-0j.- ^yj^ ) 


■'dr 


^dr 




where Zdr is the differential reflectivity in linear units 


. ^DR ~ 


So the error factor bound of the measured reflectivity, 


, IS 


1±2a/7 


(43) 


1 + 


Pm 

J zdr 


and similarly for the vertically polarized reflectivity: 


(44) 



1 + 2v/7 cos('F^^ - 20yy -pT+dyH )+ 2p^jy ^Ixz^, cos('F^^ - 26yy 
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Figure 28. Estimated reflectivity error bounds as a function of antenna polarization isolation. Solid 
lines are for rain ( = 0.98) and dashed lines for mixed phase ( = 0.9). 


The differential reflectivity, Zdr^ is the ratio of the H and V reflectivity factor (assumed that 
to simplify): 


^HV ^VH 


1 + 2a/ 7 cos + 9 j)+ 2 Pijy 




"dr 


"dr ^dr 


1 + 2y[l 26yy 9j^ + 9yfj ]+ 2p^y X 29yy + 9 jjy ) 


- ■ dp ' ^VH ) ' ^dr'^^'^y^dp 

and so the Zdr error factor range without phase control over the transmitted field, 


(47) 
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Figure 29. Estimated differential reflectivity error bounds as a function of antenna polarization 
isolation for rain (solid; = 0.98) and mixed phase precipitation (dashed; = 0.8). 


0 =6 =6 =0 

Note that if the antenna leakage and propagation phase shifts are the same ( hv w ^ 

if the differential phase shift is small, then the error can he dramatically reduced hy transmitting the V and 
H field components 90 deg out of phase: 


(48) 
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Differential 



and Specific Differential Phase 



are defined as: 
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and 


( 51 ) 

where r is range. 




dr 


The measured Differential Phase can he expressed as a sum of the measurement error and 




'F 


dp . 


+ AW 

( 52 ) dp ^ ^ 

Similarly, the measured can he written as: 
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By eliminating the second order terms and terms without and replacing 


LDR 


SHH^Sy^) = ^H^i^r- , where N is the number of independent samples averaged and LDR is the 


N 


linear depolarization ratio, the differential phase expression can he simplified as: 
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If we can assume that the antenna cross-pol and co-pol phase shifts are the same, such that 

~ 1 and « 1 then 
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Assuming that 2^1 << 1 and 2 . << A , we can further simplify as 
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This result indicates that the Specific Differential Phase '' ^ error hounds are primarily affected hy the 

target linear depolarization ratio (LDR) and averaging (N) as long as the antenna isolation is fairly good 

{241 «l). 


\xy 1 

In most conditions i will he small such that COS * i , so the error is most likely positive: 

\ 4. 
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Figure 30. Estimated Kdp error bounds as a function of antenna isolation, using Equations (53 and 

(58. An LDR or 25 dB was used in the calculations. 




Antenna Isolation (dB) 


Eigure 31. Estimated attenuation error bounds as a function of antenna isolation, based on the 
results of Eigure 30, using Kdp (Equation (28 - solid lines) and Z-Zdr based estimators (Equation 

(29 - dashed lines). 
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The conclusion of the above analysis is that polarimetric antennas for airhome weather radars should have 
at least 25 dB (30 dB preferred) polarization isolation to he useful for estimating attenuation using the 
various retrieval methods. Kdp-hased estimators can he effective with less isolation. The V and H 
components of the transmitted field should he phased close to 90° for best measurement quality. 

4 Multi-Frequency Radar (MFR) 

The task of adding another frequency to the current X-band AWR channel is significantly more complex 
and potentially more expensive than adding radiometric or polarimetric measurement capability. A 
second radar channel will likely require a separate antenna and radar electronics, although a dual 
frequency antenna, specific for airborne use, might be easier to develop than a lightweight, flat -plate dual- 
polarized X-band antenna. But most likely, only the data system and power supply sections could be 
shared. So why consider a second frequency? 

The two most important frequency dependent factors influencing weather radar measurements are that 
backscattering efficiency of small scatterers increases as f (Equation (80)) and attenuation coefficient as 
f (Figure 48). Attenuation becomes more important with range so long-range weather radars are 
predominately lower frequency (relative to X-band) S-band (10 cm wavelength) in the US and C-band (5 
cm) in Europe. On the other end of the weather radar spectrum, cloud research radars are usually Ka (1 
cm) or W-band (3 mm) to maximize sensitivity at close range. A multi-frequency weather radar has two 
primary benefits: 

1) A greater variety of weather conditions could be effectively measured. When the total path 
attenuation is high (due to long range and/or high attenuation rate), the lower frequency radar channel will 
be more effective, while higher frequency radar should have better sensitivity to small scatterers, when 
attenuation is less of a factor. Of course the availability of quality radar components and transmitter 
power are also important factors determining sensitivity at various frequencies. 

2) The difference between radar parameters measured at two (or more) frequencies can be used to 
estimate cloud and precipitation properties. Particularly differential attenuation, derived from differential 
reflectivity gradient, can be highly correlated with attenuation and liquid water content, and can be used 
to better characterize weather hazards. 

This section investigates the potential use of multiple frequencies with AWRs. First, sensitivity, back 
scattering efficiency and attenuation is investigated at a function of frequency. Then practical system 
parameters are considered along with a model of various cloud, drizzle and rain conditions, to simulate 
the frequency dependence of radar sensitivity, attenuation and measurement range in various weather 
conditions. 

MFR Measurement Concept 
4.1.1 Point Scattering 

The transmitted power density, S, of a radar can be expressed as 


(63) ' ^(W/m")^ 

p 

where ' is the transmitted power in (W), G is the antenna gain and ^ is the distance to the point target 
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reflector in meters. 


The power density captured by a point target and re -radiated towards the radar is 


5 - 


(64) 4;zr" {Anfr^ ,(W/m")^ 


The received power at the antenna feed is the product of the antenna effective aperture (A^^) and the 
incident power density of the scattered field according to 


P = AS=A 


( 65 ) 

and since 




( 66 ) 


A^j^=GA^/47t 


the received power can be written in the standard form of the radar range equation for a point target as 


( 67 ) 


' ( 47 rfr^ 


,(W). 


The actual antenna aperture area, A, can be expressed as 
so the radar range equation can be rewritten as: 


A — A fffl i , n 
e# where 


is the aperture efficiency 


( 68 ) 


P = 


P,A^n\c 


4.1.2 Volume Scattering 

The radar range equation for a point target can be extended to a volume target by replacing ^ with the 

product of volume reflectivity, / (tn / m ^ and the volume of a range cell, V. The range cell volume of 
a pulsed radar, with an approximately Gaussian shaped antenna pattern, may be approximated as 

71 7.71 

y = r" — J J f{e,^)%m{e)ded^ = r" /81n2 

( 69 ) ^00 2 

where, 

r is the range to the center of the volume cell in meters, ^ and ^ are azimuth and elevation angles 
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0 d 

respectively, and are the half power antenna heamwidths in azimuth and elevation, 


(70) 


^ , 41n(2)^\ ^ 41n(2)^% 

f = ^2 )e^( ^ ) 


e. 


3dB 


(f). 


MB 


c — 3 X 10 jn/s, speed of light in free space, 

^ “ transmit pulse length (sec.). 

For most weather radars, the antenna half power heamwidths in elevation and azimuth may he 

0 = d) 

approximated to he equal so (radians). 

The resulting expression for the received power from a volume target is: 
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This result shows that the radar received power from a volume target is a function of the power-aperture 

ct/ 

product, range squared, range resolution ( and volume reflectivity of the target. Volume reflectivity 

is dependent on frequency, polarization and beam pointing direction. This hackscattering efficiency of 
small scatterers is dramatically higher (f) at higher frequency, as shown next. Note that the Signal to 
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Noise Ratio, SNR, usually increases as range resolution squared, because the noise power decreases with 
increasing pulse length (decreasing bandwidth) and thus range resolution. 

The Volume reflectivity can be written as a sum of all the individual point scatterers in a unit volume: 


(77) 



If each scatterers can be approximated as a sphere, then the radar cross section can be estimated as the 
physical cross section 4 and an efficiency coefficient: 
a = ^ 

(78) 4 

where ^ is the back scattering efficiency of a single scatterer and D is the diameter. 



Figure 32. Radar backscattering efficiency ( ^ ) of a metal sphere as a function of radius (a) and 

radar wavelength ( ^ ) (from: [3]). 
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In the Rayleigh region, 


^-4 


(79) 

so 
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where ^ + 1 , and n is the complex index of refraction. 


at 9.4 GHz 


K\ ~ 0.97 


; at 35 GHz 


/«: !^ 0.88 


and at 95 GHz 




In the optical region, 4 ~ ^ so 4 ,whenbl»l. 

In the Rayleigh region: 
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In the Optical region: 
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4.1.3 Attenuation 


The attenuation ( ^ ) of liquid clouds and precipitation (liquid Rayleigh scatterers) is inversely 
proportional to wavelength squared and proportional to liquid volume [24] : 


Liquid in the Rayleigh region: 


(83) 


a cc 





Attenuation of ice clouds is much lower and weaker related to wavelength: 


Ice in the Rayleigh region: 


(84) 
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4.1.4 MRF Measurement Parameters 

Ideally we would like to know the particle size distribution f(D) 

, hut radar can only measure integrated 

parameters: 

<x \D^f(D) 

Reflectivity in the Rayleigh region: J -> severely biased by the largest particles 

\D^f{D) 

Reflectivity in the Optical region: J -> areal cross section 

oc r v(l>) 

Attenuation in the Rayleigh region: •> -> liquid volume 
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Number density (Nv) is also an important physical parameter. 


Let’s assume a narrow size distribution such that 


^ I 'A',, ^ 1^,12..-. 


' ^ui=\ i n4 \ \ ^Ui=\ ‘ n4 \ \ V ( 

( 85 ) " ' A 

for Rayleigh scatterers, 
and 


Z A',, ^ r^2 2 

O' ~ — > D =—N,,D 
( 86 ) 4 4 

for Optical scatterers. 

where is an equivalent particle diameter. 
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4.1.5 MFR Information Content 


Two frequencies - both in Rayleigh and no differential attenuation (high altitude ice and ash with two 
radars): 


( 87 ) 
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Two linearly dependent equations, so there is no additional information by having two frequencies: 


( 89 ) 



No New Information! 


Two frequencies - one Rayleigh, one Optical and no differential attenuation (high altitude ice and ash 
with radar and lidar sensor): 
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( 90 ) 
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Two independent equations, two unknowns {D^ and Ny). 


(92) 



Size Estimate 


From either power measurements, Ny can be also estimated. The above two cases (Rayleigh-Rayleigh 
without attenuation and Rayleigh-Optical) is well illustrated in Figure 41. The small, at both frequency 
(Rayleigh), cloud droplets form a flat line at 0 dB F-Zdr for the full range of LWC indicating that those 
radar measurements have no information content regarding LWC. On the other hand, the drizzle and rain 
points form an arcing line with high degree of correlation with liquid. Only scatterers in the Rayleigh to 
Optical transition, called Mie region would be difficult to measure. 


Two frequencies - both in Rayleigh with differential attenuation (super-cooled liquid clouds and drizzle): 
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(96) ^ ^ ^ ^ — > Liquid Water Content Information 

From either power measurement, Ny and De can also be estimated. 

In practice, retrievals are more complicated because the drop size distributions can be broad, so 
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reflectivity measurements are dominated by larger particles than attenuation in the Rayleigh region. Also 
higher frequency radars will often operate in the Mie region (Rayleigh to Optical transition), where 

with wider drop size distribution, an approximate measure of drop size and number density can be 
estimated and used to improve hazard detection and characterization. Clearly, there is potential for 
additional information from dual frequency radar when the differential attenuation rate is sufficient to 
decouple the measurements. The logical next questions are: how much frequency separation is needed; 
and what is the overlapping detection range (sensitivity) of various frequency combinations, based on 
currently available radar components? These questions are addressed next. 

MFR Components - Current Technology 

4.1.6 Antennas 

Airborne weather radar antennas are almost exclusively flat plate slotted arrays as shown in Figure 33. 
These antennas combine exceptional mechanical and electrical properties. Their light weight and narrow 
profile make them ideal for airborne use, the waveguide construction minimizes loss and thus improves 
efficiency and the absence of feed blockage leads to very low sidelobes and high beam efficiency. In size 
they range from 12” on the smaller commuter aircrafts to 20” and 30” in diameter on the mid- to larger 
aircrafts. This antenna type has been manufactured up to Ka-band frequency, but only available in single 
polarization. Dual-polarized microstrip arrays have been designed but their limited power handling 
capability makes them less attractive with weather radars. 



Figure 33. X-band flat plate slotted antenna array. 

For the following multi-frequency radar simulations we will consider the 20” diameter antenna at the 
frequencies shown in Table 1. 
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Table 3. Approximate antenna half-power beamwidth and gain as a function of antenna diameter 
and frequency. Further analysis of this report will consider the optimistic 20” diameter antenna 

size for all frequencies. 


Diameter 

X 

Ku 

Ka 

W 

12” 

7.5° (27 dB) 

4.7° (31 dB) 

2.0° (38 dB) 

0.75° (47 dB) 

20” 

4.5° (31 dB) 

2.8° (35 dB) 

1.2° (43 dB) 

0.44° (51 dB) 

30” 

3.0° (35 dB) 

1.9° (38 dB) 

0.8° (46 dB) 

0.30° (55 dB) 


4.1.7 Transmitters 

Airborne weather radars primarily use magnetron power sources for their small size, high peak power. X- 
band magnetrons in the 10-25 kW peak power range have been manufactured in high volume for marine 
ship detection and airborne weather radar applications, resulting in per unit cost below $1000 - by far the 
least expensive high power source per unit output power. Figure 34 illustrates the peak power of various 
power sources as a function of frequency. Magnetrons are vacuum devices, available up to Ka-band at 
peak power ranging from about a IkW to hundreds of kW, and 0.1% maximum duty. Other, usually 
higher average power vacuum tubes include Klystrons, Gyrotrons and Traveling Wave Tube amplifiers, 
but these typically require much larger modulators and power supplies, making them less attractive for 
airborne use. High duty solid state devices and microwave modules are sufficiently compact for airborne 
use and over the past decade have become available with average power comparable or exceeding that of 
magnetrons, but their peak power is low, and it is difficult to utilize more than 5% duty cycle with a 
single antenna radar system. Table 4 through Table 7 lists the various currently available transmit sources 
at various frequencies. For the sensitivity calculation 10 kW peak power was assumed at X, Ku and Ka 
bands and the optimistic (due to cost, size and weight) 1 kW Klystron peak power at W-band. 



0.1 1 10 100 1000 


Frequency (GHz) 


Figure 34. Transmitter power versus frequency (from the Litton-L3 MPM web site). 
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Table 4. X-band transmitters (8 - 12 GHz). 


Type 

Peak 

Duty 

Practical 


Cost 

($K) 




Peak 

Average 


Magnetron 

10-1000 kW 

0.1% 

10-20 kW 

10-20 W 

1 

TWTA 

0.1-150 kW 

CW - 6% 

- 

- 

100-1- 

Klystron 

0.5 - 3 kW 

CW - 5% 

- 

- 

100-1- 

MPM 

10-100 W 

5% - CW 

100 w 

low 

10 

Solid State PA 

20 W 

5%-CW 

20 W 

2 W 

10 

Solid State Array (20x20) 

0.1 W/element 

5%-CW 

40 W 

4 W 

10 


Table 5. Ku-band transmitters (12 - 18 GHz). 


Type 

Peak 

Duty 

Practical 


Cost 

($K) 




Peak 

Average 


Magnetron 

10-1000 kW 

0.1% 

10-20 kW 

10-20 W 

1 

TWTA 

0.1-60 kW 

0.1-30% 

- 

- 

100-1- 

Klystron 

0.5 - 3 kW 

5% - CW 

- 

- 

100-1- 

MPM 

10-100 W 

5% - CW 

100 w 

low 

10 

Solid State PA 

20 W 

5% - CW 

20 W 

2 W 

10 

Solid State Array (20x20) 

0.1 W/element 

5% - CW 

40 W 

4 W 

10 
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Table 6. Ka-band transmitters (27 - 40 GHz). 


Type 

Peak 

Duty 

Practical 

Cost 

($K) 




Peak 

Average 


Magnetron 

10-100 kW 

0.1% 

20 kW 

20 W 

50 

TWTA 

0.1-60 kW 

0.1-30% 

- 

- 

100-1- 

Klystron 

2kW 

5% 

2kW 

100 w 

100-1- 

MPM 

10-100 W 

5% - CW 

100 W 

low 

10 

Solid State PA 

20 W 

5% - CW 

20 W 

2 W 

10 

Solid State Array (20x20) 

- 

- 

- 

- 

- 


Table 7. W-band transmitters (75 - 110 GHz). 


Type 

Peak 

Duty 

Practical 


Cost 

($K) 




Peak 

Average 


Magnetron 

- 

- 

- 

- 

- 

TWTA 

0.05 - 3 kW 

10 % - CW 

- 

- 

100-1- 

Klystron 

1 kW 

5% - CW 

1 kW 

50 

100-1- 

MPM 

10-100 W 

5% - CW 

100 W 

low 

10 

Solid State PA 

20 W 

5% - CW 

20 W 

2 W 

10 

Solid State Array (20x20) 

0. 1 W/element 

5% - CW 

40 W 

4 W 

10 
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MFR Detection 


A cloud and precipitation model was used to prediet deteetion (unity single pulse signal to noise ratio) 
using the following assumed radar specifications: 


Table 8. List of radar parameters used in the sensitivity calculations. 



X 

Ku 

Ka 

W 

Peak Power (kW) 

10 

10 

10 

1 

TX Loss (dB) 

0.25 

0.3 

0.5 

1 

Antenna Gain (dB) 

31 

35 

43 

51 

RX NF (dB) 

1 

1.5 

3 

6 

RX Bandwidth (MHz) 

2 

2 

2 

1 

RF Parts Cost (2013) 

10k 

7 

100k 

300k 


Transmitted peak power: 10 kW at X, Ku, Ka and 1 kW at W-band. 

Transmitter Loss: 0.25, 0.3, 0.5 dB and 1 dB at X, Ku, Ka and W-band respeetively 
Receiver Noise Figure: 1.0, 1.5, 3.0 and 6 dB at X, Ku, Ka and W-band respectively 
Range Resolution: 150 m (1 MHz pulse bandwidth) 

Receiver bandwidth: 2 MHz with X, Ku and Ka (twice the pulse bandwidth with magnetron radars); 1 
MHz with W-band Klystron transmitter. 
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Figure 35. Single pulse, unity Signal to Noise Ratio (SNR) in clear air with 20” diameter antennas. 


4.1.8 Cloud and Precipitation Model 

In this section a cloud and precipitation model is presented for the evaluation the measurement range and 
sensitivity of various radars considered in 4.1.7. The model consists of groups of three types of size 
distrihution: 1) Clouds using the Modified Gamma distribution model, 2) Drizzle using a modified 
Marshal-Palmer distrihution model and 3) Rain using the Marshall-Palmer distrihution [pp. 306 - 317, 
Ulahy et. al, 1984: Microwave Remote Sensing]. 

Cloud model - Modified Gamma drop size distrihution: 

The number of drop per cubic meter per drop radius unit interval 

( 97 ) p{r) = ar" Qyi^[-br^) 
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where 


r = drop radius, m, 

^ and ^ are shape parameters, 




Z? is a constant related to the mode radius, ^ in units of m. 


b = 


(98) 


a 

rlr 


m„ 


and a is a function of the liquid water content, '' , according to 


a - 


(99) 

( 100 ) 

and ) is the Gamma function. 




9 

« + 4 


7 


For this simulation the parameters were varied: 

5<a<7 

Drizzle and Rain model : 

The number of drop per cubic meter per unit radius interval 

(101) p{d) = N^e.y:^{-bd) 

where 

d is drop diameter in units of meters, 

^-30x10^^. ^ = 5700/?;-' 

A^.=8xl0^ m-,^ = 4100/?;“%,„am, 

D 

and rain rate, ^ was varied from 0.1 to 5 mm/hr for drizzle and 0.1 to 150 mm/hr for rain. The resulting 
range in Liquid Water Content as a function of rain rate is illustrated in Figure 36. 
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Figure 36. Drizzle and Rain data distribution of Rain-Rate vs. Liquid water Content. 

A data base of drop size distributions was generated using tbe three models for cloud, drizzle and rain. 
This simulated data set was then used to estimate the measurement range and sensitivity of various radar 
configurations. The distribution of mean diameter and the sixth moment (mean Z diameter) of the 
distributions is shown in Figure 37. The deviation of the scatter plot from the straight x=y line is an 
indication of the spread of the drop size distribution and how the radar reflectivity is dominated by the 
largest drops in a volume cell. The data set distribution of radar reflectivity vs. attenuation at C, X and 
Ka-band is shown in Figure 38, Figure 39 and Figure 40. 
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Figure 37. Combined Cloud, Drizzle and Rain data distribution of Mean Z diameter vs. Mean 
Diamater - an indication of the dominant drop diameter of radar reflectivity factor. 
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Figure 38. Combined Cloud, Drizzle and Rain data distribution of radar reflectivity factor, Z, vs. 

Attenuation at C-band (5 GHz). 
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X — band 



Attenuation (dB/km) 


Figure 39. Combined Cloud, Drizzle and Rain data distribution of radar reflectivity factor, Z, vs. 

Attenuation at X-band (9.4 GHz). 
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Ka — band 



Figure 40. Combined Cloud, Drizzle and Rain data distribution of radar reflectivity factor, Z, vs. 

Attenuation at Ka-band (35 GHz). 
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The potential information content of differential frequency reflectivity measurements is illustrated in Figure 41. The 
cloud particles are Rayleigh scatterers at both frequencies so their differential reflectivity factor is 0 dB, irregardless 
of LWC - forming a line at 0 dB and indicating that LWC cannot be estimated; as was predicted in Equation (89. In 
the Mie (Rayleigh to Optical transition) region at Ka-band, resonnance can actually increase the radar reflectivity 
factor at Ka band compared to that at C-band (negative F-Zdr), but LWC (and/or drop size) estimation is still 
difficult because the same differential reflectivity can indicate multiple solutions to LWC. But when the higher 
frequency is in the optical region (the arc of points from [0 g/m^, -0.75 dB] to [3 g/m^, 8 dB]), LWC can be easily 
estimated, confirming (92. 


C — Ka Differential Reflectivity 



Liquid Water Content (g/m^) 


Figure 41. Combined Cloud, Drizzle and Rain data distribution of frequency differential 
reflectivity factor, F-Zdr, vs. Liquid Water Content (g/m^). 
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4.1.9 Detection to 10 and 20 km Range 


In this section, detection of the various simulated cloud, drizzle and rain cases are determined for the 
radar configurations (frequencies and system parameters) tabulated in Table 8 and assuming a 150 m 
range resolution. The attenuation rate and reflectivity was calculated for each cloud, drizzle and raindrop 
size distribution. Then the radar received signal to noise ratio was calculated for each simulated weather 
condition, considering reflectivity and the two-way attenuated propagation to 10 km and 20 km (the 
weather condition, and the corresponding attenuation rate, was assumed to be uniform over the two-way 
propagation distance). Detection was assumed for unity or greater signal to noise ratio. This detection 
threshold is reasonable because averaging was not accounted for, which reduces the false alarm rate. 
Undetected weather events (simulated data base points) are indicated by small black points through 
Figure 42 to Figure 44. 

Multi-frequency detection to 10 km range is shown in Figure 42 and Figure 43, and to 20 km in Figure 
44. The scatter plots clearly show that the lower frequency radars are best suited for the heavier 
precipitation, while the higher frequency millimeter wave systems have better sensitivity at short ranges. 
But even from a range of 10 km, attenuation severely diminishes detection at Ka and W-bands. As 
expected, X-band is the best compromise between sensitivity to small scatterers and attenuation, but C- 
band is the only of the tested frequencies that can penetrate the heaviest precipitation cases. 

The results indicate that detection of the most hazardous conditions like high LWC cloud and severe rain 
is not reliable at multiple frequencies. The combination of C and X bands (red and orange points in Figure 
42 and Figure 44) is effective in most drizzle and weak to moderate rain events, but C-band lacks 
sensitivity in most clouds and small particle drizzle and X-band attenuation is too strong in the heavier 
rain cases. The X and Ka combination (red and purple points in Figure 42 and Figure 44) can detect more 
cloud cases and the weaker drizzle, but ineffective in anything stronger than moderate rain. Also, from a 
range of 20 km, the attenuation at Ka-band severely reduces the number of detected cases. 

A better reason for multi-frequency is the option to switch to a lower frequency channel when attenuation 
at X-band prevents detection. This improvement in sensitivity at lower frequency in highly attenuation 
conditions is well illustrated in Figure 45, Figure 46 and Figure 47. C-band can detect all but the weakest 
drizzle and rain conditions to 20 km but very few of the clouds. X-band detects more of the clouds and 
just about all the drizzle but is unable to penetrate heavier rain. Attenuation at Ka-band severely limits 
detection to just a few low LWC cases. 
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MFR Detection at 10 km 
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Figure 42. C, X, and Ka-band simultaneous, multi-frequency detection to 10 km range. 
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MFR Detection at 10 km 
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Figure 43. Simultaneous MFR detection to 20 km range at X, Ka, and W-bands. 
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Figure 44. Simultaneous MFR detection at C, X and Ka-bands, to 20 km. 
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C Detection at 20 km 



Figure 45. C-band detection to 20 km range. 
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Figure 46. X-band detection to 20 km range. 
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■ Detection 
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Figure 47. Ka-band detection to 20 km range. 


MFR Differential Attenuation 

In severe weather conditions X-band attenuation can be quite high, resulting in underestimating 
dangerous weather cores, but worse: can lead to the received signal diminishing below the detectable 
threshold. So by shifting to a lower frequency, attenuation quickly decreases. The attenuation rate (dB 
km ' mm ' hr) in rain as a function of radar frequency is illustrated in Figure 48. Attenuation is 

approximately proportional to ' , so a small change in operating frequency, that is available within a 
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waveguide band, yields very small change in attenuation rate and thus is not worth the added cost and 
system complexity. However, if a C-band channel was available in heavy (50 mm/hr) precipitation, the 
attenuation rate decreases from 1.5 dB/km at X-band to 0.3 dB/km. Over a 10 km wide precipitation (20 
km two way propagation) this adds up to a 24 dB difference in total attenuation, and approximately a 12 
dB improvement in sensitivity; assuming similar radar power-aperture product. So when X-band 
attenuation is so severe that attenuation correction techniques, such as polarimetric or radiometric, 
cannot be employed, C-band radar might still penetrate heavy precipitation and detect and identify 
hazardous weather cores. 



Figure 48. Attenuation rate per rain rate as a function of radar frequency (from [3]). 


There is a potential benefit in adding a second, higher frequency. The backscattering efficiency of small 

f4 

(Rayleigh) hydrometeors increases as ' , so for example a Ka-band (33 GHz -> 3.5 times that of a 9.4 

GHz X-band) radar, given the same radar power-aperture product, receiver sensitivity and attenuation, 
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will have almost 22 dB more sensitivity. The problem however is attenuation. Attenuation can easily 
overcome this gain. The attenuation rate at X-hand is approximately 0.006 dB km ' mm ' hr while at Ka- 
hand is 0.3 dB km ' mm"' hr. So even a 10 mm/hr rain will cancel the hackscattering sensitivity gain over 
a short, two-way propagation distance of 7.5 km (3.75 km range)! Shifting from X-hand to Ka-hand only 
makes sense when attenuation is low; i.e. for the detection and mapping of freezing drizzle/light -rain or 
high altitude ice clouds. 

For the same reason, dual-frequency X-Ka hand measurements also have limited use. Although the 
differential attenuation rate is sufficiently high for easy detection and estimation of rain rate and 
attenuation, as shown in Figure 52, in those highly attenuating conditions, dual-frequency measurement 
range is severely limited as confirmed hy the results shown in Figure 42, Figure 46 and Figure 47. 

The following calculation finds the measurement accuracy that is necessary for useful dual-frequency 
measurements: 


( 102 ) 

(103) 

and 

(104) 


(105) 


P =P +P 

m s n 

P +P 

cr,„ = 


Jn . 
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The standard error deviation of a reflectivity measurement in units of dB can he expressed as 


dBE^ =l0log 


10 


V P. y 


= 101og 


10 


(106) 

In high SNR this expression simplifies to 
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The estimation of differential attenuation ( ) requires four reflectivity measurements: 


M/,./2) = 


{dBZ,^ fc ) - dBZ,^ (r, ))- {dBZj^ (r, ) - dBZ,^ (r, )) 


(108) 


r^-r, 


dBE 

Consequently, the expected error in the measured differential attenuation, ^ , in high signal to noise 
ratio conditions can he expressed as 


(109) 


dBE^ 


x /4 X dBE ^ 
Ar 


If we are satisfied with a 1 km range resolution and a 1 dB/km differential attenuation as a practical goal, 
then the measurement standard deviation should he no more than 0.5 dB. To achieve this accuracy, the 
required number of independent samples is 


( 110 ) 


0.5 = 101og 


10 


1± 


1 

4n , 




Note that the resulting measurement quality is still poor since the 0.5 dB error corresponds to a single ^ 
deviation, yielding only a 68% of the measurements within 0.5 dB of the actual differential attenuation, as 
shown in Figure 49. To increase the confidence interval to 95%, 287 samples are required. Another way 
to look at this is that if the differential attenuation is 1 dB, then our measurements should have much 
smaller than 1 dB standard deviation. If 287 samples are averaged, then the resulting differential 
attenuation measurement standard error (deviation) will he approximately ±0.25 dB. But if -300 
independent samples are required for a single differential attenuation measurement, then a complete 
volume scan ahead of an aircraft would take a long time. 


If an AWR has a PRF of 5 kHz (30 km maximum range), then approximately 17 (5000/300) 
measurements can he made in a second, so a 100 point volume scan would take about 6 seconds. But this 
calculation assumes that all the samples are independent at all the frequencies. At X-band the approximate 
decorrelation time of precipitation due to differential motion of the particles, is -10 ms. A rate 50 times 
slower than the 5 kHz PRF assumed above. So either much longer integration time is required, or each 
radar has to be equipped with frequency agility to force decorrelation by shifting the operating frequency 
by at least the transmitted pulse bandwidth [25] . Even if all these obstacles are overcome, what frequency 
combinations and what weather conditions would have 1 dB/km or greater differential attenuation? 
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Normal Curve 
Standard Deviation 



Figure 49. Gaussian distribution as a function of standard deviation, a. The ± (J interval 
corresponds to about 68% of the samples, while +2 (T corresponds to 95% and ±3 a to 99.7%. 


The cloud, drizzle and rain model established in Section 4.1.8 is used again to investigate differential 
attenuation of various frequency combinations. The results are shown in Figure 50 through Figure 52. 

Only the heavier rain produces a differential attenuation rate more than 1 dB/km at X-C and X-Ku 
frequency combinations. The X-Ka frequency combination is sufficient in most rain and higher LWC 
cloud conditions, but attenuation also diminishes the sensitivity of a Ka-band radar at longer ranges as 
shown in Figure 50 and Figure 5 1 . 

We can conclude that useful differential reflectivity measurements are quite difficult. The differential 
attenuation needs to be quite high, often resulting in impractically short maximum detection range at the 
attenuating frequency, a large number of independent samples are required, and the resulting retrievals 
have coarse range resolution. 


73 



2 - 


k: 

< 

E 

o 


1 - 


0 
0.1 


X - 


C Differential Attenuation (dB/km) 

-TTTl- 


T 


T 


3 dB/km 






. »»«- >»#» 

A --t* »<•--*** 


V. 

' * 

4.. ' * ^ |c ' i 

^w» 

* - 5- _ 

-;.ivc?tfy::, 

*jT-j» * ^ ■ f#+-? 

♦ T , . *r* 

’'4^'” n* i\ ,, 

▼ -- 


t , *- 


lL 


1.0 10.0 100.0 

Mean Diameter (/:im) 


I rr i 



1000.0 


Figure 50. X vs. C-band differential attenuation for a range of cloud, drizzle and rain conditions. 
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Figure 51. Ku-X differential attenuation estimates are less than 0.5 dB/km for clouds and drizzle, 
and only reach 2 dB/km in the heavier rain conditions. More frequency separation is needed for 

MFR retrieval. 
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Figure 52. Ka-X Differential attenuation. 


5 Sensitivity 

A limitation of X-band weather radar is often the lack of sensitivity due to either low scatterer cross 
section (e.g.: high altitude ice crystals or volcanic ash), high attenuation (severe thunderstorms) or a 
combination of both (freezing drizzle and rain). The only remedy to attenuation is to shift to a lower 
frequency, but improving sensitivity through adjusting the radar operating parameters and improved 
signal processing will always improve the probability of detection. A few of these methods are presented 
next. 

Following is the explanation why the easiest and quickest way to boost the sensitivity of a weather radar 
is to degrade its range resolution. Low radar cross section is primarily due to the small size of scatterers, 
but low index of refraction can also be an important factor (glass composite of volcanic ash). From (76), 
the received power from a volume scatterer can be expressed as: 
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p. = 


n 


( 111 ) 




32r" 


The corresponding thermal noise floor, Pn, is: 


P„ - kTBF = 

( 112 ) ^ , 

where k is Boltzman’s constant (1.38E-20 mW/K*Hz), T is temperature is deg. K, F is the receiver noise 
figure and B is the receiver bandwidth (in units of Hz) matched to the transmit pulse length ^ . 

Signal to Noise Ratio (SNR) is a key factor of detection, in addition to averaging and setting a threshold 
level to an acceptable probability of false alarm rate. SNR is a ratio of the measured signal and receiver 
thermal noise floor: 


PAeifCT 

cyn _ Pr 32r^ .. 

P„ kTF 32r^kTF 

( 113 ) ^ 

Note the square of the transmit pulse length in the numerator of the expression for SNR. This indicates 
that for a given weather radar with fixed peak transmit power, antenna aperture and noise figure, 
degrading range resolution, i.e. transmitting longer pulses, is the most effective way to improve 
sensitivity. With a receiver bandwidth matched to the transmitted pulse, the measurement signal to noise 
ratio improves as square of the range resolution. 

Another method of improving sensitivity is through signal processing. In low turbulence regions, a 
combination coherent (FFT) and non-coherent (power) integration can boost the probability of detection 
without lowering the false alarm rate. The comparison of conventional non-coherent averaging and a 
combination of coherent (FFT) and non-coherent averaging is illustrated in Figure 53. A W-band ground- 
based radar was used to collect data from a non-precipitating cloud on February 9, 2012 in Amherst MA. 
The radar was configured to transmit at constant PRF and data blocks from 10240 transmit pulses were 
processed three ways: 1) the top figure shows the combination of FFT and non-coherent averaging, where 
the 10240 data points where broken up into groups of 128 for FFT processing and the 80 power 
spectrums where averaged (non-coherently); 2) the middle figure illustrates the conventional power 
averaging of all 10240 data points and the 3) bottom figure is commonly called Coherent Power, 
generated by forming 620 single lag pulse -pairs and averaging those complex correlations and displaying 
the magnitude of the resulting average. 

The false alarm rates are quite similar in all three methods, but the FFT shows the best detection, followed 
by the power average and the coherent-power close third. The Coherent Power image could have been 
improved by forming 10239 pulse pairs from the 10240 pulse train, which likely would have resulted in 
an image quality similar to that of the power average, but with more uniform false alarm background. 

An even better method of detecting weak scatterers is using the Doppler velocity image. Figure 54 shows 
the same data but processed to velocity using the FFT processing (top) and pulse-pair (bottom). The 
coherent integration of the FFT processing still has the best detection, but even the pulse-pair velocity 
image surpasses all three power image detection methods of Figure 53. 
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Figure 53. Top: Feb 9 09:00 UTC W-band co-pol FFT-based dBZ (top); non-coherent dBZ 
(middle); pulse-pair lag 1 dBZ (bottom). Conditions: pfa=0.5; PRF: 4.960 kHz; FFT: N=128; 
M=80; non-coherent and pulse-pair lag 1: NM=10240. Data gathered at an elevation angle of 145 

degrees. 
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Figure 54. Top: Feb. 9 09:00 UTC W-band co-poI data FFT-based velocity (top); PP velocity 

(bottom). 
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